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ABSTRACT 

We revisit the notion that galaxy motions can efficiently heat intergalactic gas in 
the central regions of clusters through dynamical friction. For plausible values of the 
galaxy mass-to-light ratio, the heating rate is comparable to the cooling rate due to 
X-ray emission. Heating occurs only for supersonic galaxy motions, so the mechanism 
is self-regulating: it becomes efficient only when the gas sound speed is smaller than 
the galaxy velocity dispersion. We illustrate with the Perseus cluster, assuming a 
stellar mass-to-light ratio for galaxies in the very central region with the dark-matter 
contribution becoming comparable to this at some radius . For < 400 kpc ~ 
Srcooi — corresponding to an average mass-to-light ratio of ^ 10 inside that radius — 
the dynamical-friction coupling is strong enough to provide the required rate of gas 
heating. The measured sound speed is smaller than the galaxy velocity dispersion, 
as required by this mechanism. With this smaller gas temperature and the observed 
distribution of galaxies and gas, the energy reservoir in galactic motions is sufficient to 
sustain the required heating rate for the lifetime of the cluster. The galaxies also lose a 
smaller amount of energy through dynamical friction to the dark matter implying that 
non-cooling- flow clusters should have flat-cored dark-matter density distributions. 

Key virords: galaxies: clusters, general - galaxies: evolution - galaxies: formation - 
galaxies: interactions - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Galcixy cluster gas loses thermal energy copiously through 
X-ray emission. In the absence of energy input, radiative 
cooling in cores of clusters should result in substantial gas 
inflow (see Fabian 1994 for a review). These "cooling flows" 
would have associated mass deposition rates of several hun- 
dred M0yr~^ in some clusters (Peres et al. 1998). Never- 
theless, recent high-resolution X-ray observations (e.g., Pe- 
terson et al. 2001, 2003) have revealed that there is little 
evidence for the expected multi-phase gaseous structures, 
strongly suggesting that mass dropout is being prevented 
by some source that is heating the gas, thereby balancing 
radiative energy loss in the central region of clusters. Re- 
cent work has focused on two prospective heating mecha- 
nisms: (1) diffusive heat transport, via thermal conduction 
(Tucker & Rosner 1983; Bregman & David 1988; Narayan & 
Medvedev 2001; Voigt et al. 2002; Fabian, Voigt, & Morris 
2002; Zakamska & Narayan 2003; Kim & Narayan 2003a) 
and/or turbulent mixing (Cho et al. 2003; Kim & Narayan 



2003b; Voigt & Fabian 2004), from hotter gas in the outer 
region to that in the core; (2) energy input from jets, out- 
flows, and radiation from a central active galactic nucleus 
(Ciotti & Ostriker 2001; Churazov et al. 2002; Briiggen & 
Kaiser 2002; Kaiser & Binney 2003). 

There is however at least one additional mechanism that 
appears to have been overlooked. It involves the energy lost 
by concentrated clumps of matter (galaxies) as they move 
through the cluster. Part of this energy may go into re- 
arranging the dark-matter mass distribution (El-Zant et al. 
2003), but a significant fraction should end up deposited in 
the gas. That dynamical-friction (DF) coupling can trans- 
form the dynamical energy of galaxies into thermal motion 
in the gas has been known for at least four decades (e.g., 
Dokuchaev 1964). Relatively recent work involving this no- 
tion includes the analysis by Miller (1986) of the Perseus 
cluster and that of Just et al. (1990) concerning the Coma 
cluster. Both studies confirm that, provided that the mass- 
to-light ratio of galaxies is ~ 20, energy lost by galaxies to 
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the gas should be sufficient in counteracting the radiative 
coohng of the gas in the central region of these clusters. 

Several developments on the empirical side suggest re- 
newed relevancy for this mechanism. One obvious one in- 
volves recent X-ray data confirming that a heating mecha- 
nism is required; whereas the consensus in the 1980's was 
against this conclusion, it now seems inescapable. The sec- 
ond involves revisions to the inferred gas electron densi- 
ties in the central region of clusters; best values referred 
to by Miller and Just ct al. arc of the order of lO^'^cm^'', 
while current best estimates are rather in the range of 
10"^ - 10"^cm"^ (Kaastra et al. 2004). This leads to an 
ordcr-of-magnitudc increase in the dynamical friction cou- 
pling between the galaxies and gas. There has also been 
progress in determining the mass-to-light ratio of galaxies. 
On the theoretical side, work by Just et al. (1990) and Os- 
triker (1999) has since shed some light on the behavior of 
the dynamical-friction coupling in a gaseous medium in the 
transonic and subsonic regimes. 

There appears to bo a priori no reason why the rate 
of energy loss from gala^xies to gas via dynamical friction 
should be of the same order of that radiated by the gas. 
However, this coupling is active only when the sound speed 
of the gas is smaller than the typical velocity of galaxies. The 
mechanism is therefore self-regulating; the gas cools until 
the dynamical-friction heating rate is always equal to the 
cooling rate. We start by pointing out, in the next Section, 
why this is expected to be so, before moving on to develop a 
Monte-Carlo model to estimate the total energy expected to 
flow into the cooling region of the Perseus cluster, using re- 
cent data for galaxy luminosities and projected positions, as 
well as for the gas parameters of that cluster, and averaging 
over a set of different realizations where three-dimensional 
positions, velocities, and mass-to-light ratios are treated as 
stochastic quantities. The final Section discusses briefly the 
central dark-matter distribution and energetics, as well as 
some remaining questions. 



2 DYNAMICAL FRICTION IN GASEOUS 
SYSTEMS 

The fundamental difference between the dynamics of coUi- 
sionless and gaseous systems is the existence in the latter 
case of pressure gradients which communicate forces at the 
sound speed Cs ■ If the galaxy velocity V is much larger than 
Ca then the gravitational interaction between a particle at 
impact parameter b leading to dynamical friction is unaf- 
fected by pressure forces (because by the time these are com- 
mu nicated to th at region the galaxy is already at a distance 
b^/l + {V/cay > 6; e.g., Ruderman & Spiegel 1971). If, on 
the other hand, V < Cs, pressure forces can be communi- 
cated to the perturbed region before the minimum distance 
b is achieved. The resulting pressure gradients ensure that 
the displacement of gas particles therein is hindered, and so 
the cumulative back reaction on the perturber that results 
in the dynamical-friction force is drastically reduced. 

In the collisionless case, when one speaks of a parti- 
cle, one has in mind an actual material constituent. In the 
classic Chandrasekhar formulation, each of these contributes 



individually. The sum of the perturbations from particles 
moving faster than the perturber (a galaxy) is negligible for 
small particle mass (being proportional to that mass) . That 
from slower particles however results in a net contribution 
that for the massive perturber is independent of the smaller 
background-particle mass, and is always directed opposite 
to the perturber's motion. This results in the gradual de- 
crease in the dynamical-friction force when the perturber 
velocity falls below that of the typical background parti- 
cle. In the gaseous case however, unless the mean free path 
is larger than the impact parameter, collisions will lead to 
rapid decoherence of individual particle motion. A coherent 
effect that results in dynamical friction should then involve 
interactions between the perturber and gaseous elements in- 
stead of individual particles. If the bulk velocity, that is the 
macroscopic motion, of the gas elements is smaller than that 
of the perturber then all gas elements will contribute. 

The above considerations suggest that gaseous dynam- 
ical friction will follow the high-velocity approximation in a 
collisionless medium when V > Cs and then drop sharply 
when V < Cs- This conclusion was already reached, on 
the basis of steady-state perturbation theory (Ruderman 
& Spiegel 1971; Raphaeli & Salpeter 1980) and was con- 
firmed by two more sophisticated, but rather different, tech- 
niques invoking time dependence (Ostriker 1999) and fluc- 
tuation theory (Just ct al. 1990). Simulations by Sanchez- 
Salcedo & Brandenburg (1999) have shown broad agreement 
for particles moving on rectilinear trajectories and, more im- 
portantly, qualitatively confirmed the above conclusions in 
spherical systems (Sanchez-Salcedo & Brandenburg 2001). 
For a perturber of mass Mp, moving with V > Ce vfe will 
thus assume that energy is lost at a rate given by 

dE _ AiriGMpf Pg 1 &max / ^ \ 

it~ V 6~r' 

and that the energy loss vanishes for V < Cg. 

The relevant density pg is the typical gas density to be 
found within the cluster cooling radius rcooi, inside which 
the cooling time is less than a Hubble time — that is, the re- 
gion where heating must be invoked. If the galaxy is inside 
this radius, the range of impact parameters that contribute 
to heating should, to a first approximation, correspond to 
a minimal scale determined by the galaxy size and a maxi- 
mal one determined by the cooling radius. For these galaxies 
therefore we take the Coulomb logarithm to be In (rcooi/ro), 
where the "galaxy radius" ro is set to 10 kpc — about twice 
the effective radius and disk scale length for ellipticals and 
spirals, respectively. Galaxies outside TcooI should also con- 
tribute, but in this case the contribution from material inside 
''cool to the energy lost by a galaxy at radius r comes from 
impact parameters between max(ro, r — rcooi) and r -|- rcooi. 
This contribution is also restricted from within an angle 
2tan~^[rcooi/('' — fcooi)]- Replacing the minimal impact pa- 
rameter with r— rcooi leads to an error of < 1% in the results. 
For galaxies outside rcooi therefore 

In ""^^ —> — tan~ 

1 rcooi r -I- rcooi ^2) 

^min r rcooi r rcooi 

For galaxies on highly eccentric trajectories, the bulk of the 
energy exchange with gas inside the cooling radius takes 
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Figure 1. Number of galaxies inside projected radius R, which 
can be approximated by a linear fit (top) over a large range in radii 
and with limited relative error (bottom). The arrow indicates the 
cooling radius. 



place at closest approach to rcooi- However, corrections due 
to this effect are not large, even if most galaxies are in- 
deed on nearly radial orbits. For any individual realization 
of such a highly anisotropic quasi-steady-state system, the 
dominant contribution to the energy input will come from 
those galaxies that are at small radii (since both terms on the 
right hand side of ^ rapidly decrease when r is appreciably 
larger than Tcooi) and are thus on their closest approach. 



3 THE PERSEUS CLUSTER 

Because of the rather large energy input rate (~ 10*^ erg/s) 
that is required to prevent a cooling flow, the Perseus cluster 
is one of the more demanding cases for any cluster-heating 
model. In this Section we examine whether, even in this case, 
energy dissipated from cluster galaxies can be sufficient to 
provide the required power. 

Brunzendorf & Meusinger (1999) have collected a cata- 
log of projected positions and apparent luminosities of 660 
galaxies in the Perseus cluster. In Fig.^we plot the cumula- 
tive number of galaxies as a function of projected radius R. 
As Fig.0shows, this can be fitted over a large range in radius 
by a distribution N{R) ~ R, which implies a projected den- 
sity distribution n{R) ~ 1/R. This in turn implies a three 
dimensional density varying as ~ l/r^. In this case, along 
the line of sight, and at any given R, the number density of 



galaxies decreases as 



(i?2 + Z'^) 



(with r 



R^ + 



This simple dependence of the number-density distribution 
on Z is consistent with the latter being distributed in such 
a way that \Z\ = i?tan~^ X, with X being a random vari- 
able uniformly distributed between and 7r/2. This is how 
the third spatial coordinate is chosen in the random real- 
izations of the Perseus cluster that we describe below. Since 
the density does drop more sharply than at large radii, 
this prescription will slightly overestimate the Z values of 
coordinates. This is however a minor effect, since that re- 



gion does not contain many galaxies and since most of the 
DF comes from inner cluster galaxies (and if anything leads 
to a decrease in the total DF power pumped into rcooi)- 

Galactic extinction is particularly severe at the posi- 
tion of the Perseus cluster in the sky. Trials with the NED 
calculator show that corrections in the central region vary 
between ~ —0.65 and ~ —0.9, with a value of about —0.7 
magnitude at the CD galaxy NGC 1275 (which is not in- 
cluded in the DF calculations). To all galaxies in the cata- 
log we will apply an extinction correction of —0.75 magni- 
tude. To this one needs to add a K-correction that ranges 
from -0.092 for ellipticals to -0.02 for late type spirals. The 
apparent magnitude of disks are however further affected 
by inclination effects. These vary as Ai ~ — \og{a/b), with 
b/a ~ cosi (assuming the axis ratio of the galaxy viewed 
edge on to be vanishingly small) . This gives an average cor- 
rection Ai = —0.19. We also apply a cosmological dimming 
correction of 101og(l + z) = —0.085 to all galaxies. Abso- 
lute magnitudes are then obtained by assuming a distance 
to the Perseus cluster of 78 Mpc, compatible with a Hubble 
parameter Hq = 100 h Mpc with h = 0.7. 

Unless the structural parameters defining the density 
distribution of elliptical galaxies vary significantly with 
mass, the Faber- Jackson relation and the virial theorem 
imply a luminosity-dependent central mass-to-light ratio: 
M/L ~ L^, with 7 ~ 0.3 - 0.4. Gerhard et al. (2001) have 
dynamically examined the central mass-to-light ratio of a 
sample of elliptical galaxies and determined that their re- 
sults are in agreement with the assumption of structural 
homology, and that the variation of the mass-to-light ratio 
within an effective radius is largely due to change in stel- 
lar population. Although Trujillo et al. (2004) have argued 
that non-homology effects may be important in determining 
7 (which they deduce to be ~ 0.1), a stellar mass-to-light 
variation with 7 ~ 0.4 also appears to exist for disk galaxies 
(Salucci, Ashman & Persic 1991). All these relations exhibit 
significant scatter. We adopt a zero-point mass-to-light ratio 
(which we will refer to as "stellar", even though the afore- 
mentioned papers do not exclude a modest dark-matter con- 
tribution in the inner regions) consistent with these results 
and containing a stochastic component to represent the scat- 
ter. Thus for a galaxy of luminosity L we adopt 



= X 10 



5 X IQiOLp 



(3) 



where X is a random variable chosen from a normal distri- 
bution with average of unity and dispersion 0.7. (The nor- 
malization is chosen by inspection of Fig. 13 of Gerhard et 
al. 2001.) 

The mass-to-light ratio is expected to increase with 
radius, as the dark-matter contribution becomes progres- 
sively more important (e.g., Gerhard et al. 2001; Takamiya 
& Sofue 1999). For galaxies confined in clusters, the con- 
tribution of dark matter at larger radii will depend on how 
much of the halo that should be surrounding the galaxy 
has been removed. This will in turn depend on the maxi- 
mum excursion of that galaxy into the cluster center, and 
the density distributions of both the galaxy and the cluster. 
For example, a singular isothermal sphere with character- 
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Figure 2. Mass-to-light ratios for galaxies inside radius r de- 
duced from equation |3 over 10,000 random realizations of the 
Perseus cluster data. 



istic velocity dispersion a moving through another isother- 
mal system with dispersion E > cr, will be cut to a radius 
Vt = ((T/E)r,„in if rmin IS the minimum distance that the 
centers of the two spheres achieve. And so, in this approxi- 
mation of the situation, an average galaxy halo (velocity dis- 
persion ~ 200 km s~^) moving through the Perseus cluster 
(velocity dispersion ~ 2000km s~^), will have rt ~ fmin/lO 
and the mass enclosed within this radius will be (2(7^/G)rt. 
Since halos found in cosmological simulations have density 
profiles that are close to singular isothermal spheres over a 
large range in radii, this representation should be sufficient 
in approximating the initial conditions before any modifica- 
tion of these profiles takes place — provided that the density 
of the galactic halos is scaled in such a way that the mass in 
the very inner region is not dominated by the dark matter, 
so as to be in line with the aforementioned studies concern- 
ing the central mass-to-light ratios of galaxies. 

We will assume that the radial variation of the mass- 
to-light ratio of cluster galaxies is given by 

where {M/L)s is defined by equation @ and = 
GMium/cr^ ~ f(r/W) (with / = rmin/r and Mi^m the lumi- 
nous mass) defining a spatial scale over which enough dark 
matter remains tied to the galaxy for its contribution to the 
total mass to be comparable to the luminous one. Statisti- 
cally, / will be determined by the velocity distribution of 
galaxies, with more anisotropic dispersions implying deeper 
entries into the cluster core and so more stripping. Uncer- 
tainty in this parameter are represented by an additional 
appeal to another random variable Y again chosen from a 
normal distribution with mean unity and dispersion 0.7. In 
Fig.|21we show the average mass-to-light ratio within a given 
radius for the Perseus cluster averaged over 10, 000 realiza- 
tions of the random variables X and Y . 

Another dynamical variable determining the gas- 
galaxies DF coupling is the velocity distribution of the latter. 
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Figure 3. Power expected (over ten thousand random realiza- 
tions of the cluster) to be transmitted to gas inside rcool via 
dynamical friction on cluster galaxies within radius r, relative to 
the total power radiated from within rcool- 



Brunzendorf & Meusinger have only determined the line-of- 
sight speeds of 169 of the galaxies in their samples. The ve- 
locity dispersions of galaxies, when binned radially in groups 
of 20 to 30, are ~ 1000 — 1400 km s"'^ with no obvious sys- 
tematic variation incompatible with small-sample statistics. 
Fig. 13 shows the energy expected to be deposited via dy- 
namical friction into rcooi, from galaxies enclosed within a 
radius r, relative to the energy radiated from the cooling 
region, for different values of the radial scale r^. For each 
velocity component a normal distribution with zero mean 
and dispersion a, such that \/^a = 2000 km s~^, is assumed. 
Results are shown for the 10,000 realizations of the cluster 
in Fig. 131 assuming r-cooi = 130 kpc, and the total energy 
emitted from r < rcooi to be 5 x 10*'*erg/s, consistent with 
the findings of Kaastra et al. (2004) , when converted to cor- 
respond to h — 0.7. The sound speed Cs — \/ "yRT/ fj, can 
also be read off Table 5 of Kaastra et al. (2004). Excluding 
the very inner data point, which is probably influenced by 
central activity (Churazov et al. 2003), the central sound 
speed starts at about 640 km s~^, eventually rising with ra- 
dius to reach an asymptotic limit almost double that value, 
but remaining largely constant within most of the cooling 
region (Churazov et al. Fig. 8). These studies also show the 
electron density to be slowly varying within rcooi, with val- 
ues in the range 0.1 — 0.005 cm~'^ (there is also a factor of 
^y0.7/0.5 that should be included because we use h = 0.7, 
as opposed to the value h — 0.5 used by these authors). In 
calculations shown in Fig. |3] a value of 0.02 cm"'' (setting 
fj, = 1.18 to convert electron densities into gas density) and 
Cs = 700 km s~^ are used. 

Because of the sharp cutoff in the DF-mediated cou- 
pling when V < Cs, the amount of energy deposited into the 
gas must sensitively depend on the relative magnitude of the 
galaxy velocities and the sound speed when V Cs- The for- 
mer is a function of the anisotropy of the distribution. For 
example, if galaxies move on nearly radial trajectories, their 
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Figure 4. Same ordinate as in Fig.^but plotted against the gas 
sound speed for different (constant as well as radially varying) 
values of the ratio of the velocity dispersions erg / (7^ (and assuming 

velocity dispersion would be closer to the line-of-sight value 
of 1200 km s~^, rather than 2000 km s~^ as assumed above. 
The simulations of Ghigna et al. (2000) suggest that galax- 
ies (as represented by cluster subhalos) may have biased 
velocity dispersion in the central region, perhaps reflecting 
enhanced radial motion; whereas at larger radii they have 
velocity distributions that rather closely follow that of the 
dark matter. Fig. |1| shows the power pumped into the gas 
relative to that radiated within the cooling radius of the 
Perseus cluster, as predicted by our model, as a function of 
the sound speed for a range of fixed and (radially) varying 
(jfl/cTr values (assuming — ag) for Va — 400 kpc. It is 
rather remarkable that the competing effects of increased 
coupling at low velocities (because of the 1/V dependence), 
and the rapidly decreasing number of galaxies contributing 
to the DF with increasing Cs , conspire to always ensure that 
a maximum value for the energy input is approached when 
600 km s~^ <Cs < 700 km s~^, precisely the value attained 
in the (central) region where the temperature drops in the 
cluster. 



4 CONCLUDING REMARKS 

That dynamical-friction coupling between cluster galaxies 
and gas can significantly heat the latter component is not 
a new concept. Indeed Miller (1986) has pointed out that a 
few luminous galaxies in the Perseus cluster's central region 
may alone deposit sufficient energy to keep gas in that region 
from cooling, provided that they had a mass-to-light ratio 
of about 20. We revisited this cluster using recently com- 
piled galaxy and gas data and determinations of the mass- 
to-light ratios in galaxies. Our Monte Carlo model shows 
that the rate at which energy is deposited by galaxies into 
the gas within the cooling radius can completely compen- 
sate for radiative loss from within that radius, if galaxies 
at radii < 400 kpc from the center of the cluster have dark 



mass that is comparable to their stellar mass (with galax- 
ies at smaller radii having progressively smaller dark-matter 
content, vanishing for those near the center of the cluster). 
The associated average mass-to-light ratio within this radius 
is about 10. 

A robust and potentially important prediction of the 
model is that the drop in gas temperature invariably ob- 
served in the cooling regions of clusters is necessary for 
the DF coupling to be sufficiently strong. There is a nat- 
ural interpretation for this phenomenon: the gas cools un- 
til dynamical-friction heating becomes sufficiently efficient 
to prevent a further drop in temperature. The gas random 
motion is then coupled to that of the galaxies, and so re- 
mains in approximate equilibrium with it. All this simply 
follows from the fact that dynamical friction drops sharply 
for subsonic motion. In the case of the Perseus cluster, the 
transition from weak to strong coupling occurs precisely in 
the range of sound speeds found in the cooling core of the 
cluster. 

The action of the mechanism discussed in this paper 
also has consequences for the dark-matter distribution. As 
shown by El-Zant et al. (2003), dynamical friction from the 
galaxies will heat a density cusp, creating a core with ra- 
dius corresponding to roughly a fifth of the original Navarro, 
Frenk & White (1997) scale length (cf.. Fig. 1 of El-Zant et 
al.), where dark matter has been spread out to larger radii 
(Fig. 3 of El-Zant et al.). Approximating the initial density 
distribution in this region such that pi = poiro/r) and the 
final one with p/ = po, the energy required for this transfor- 
mation is A* = (SttG)-^ J^°{\d(t)/dr\'i - \d(t>/dr\j)4:7vr^dr = 
227rVo'-o/45 [with d(j}/dr = {4nG/r^) J pr^dr]. Taking ro = 
'"cool and po to correspond to four times the value de- 
fined by the electron density Pe(?"o) = 0.0053 (where we 
have used the gas to gravitational mass ratio in Table 5 of 
Peres et al. 1998 and equation (4) of Churazov et al. 2003 
for the electron density), one finds A$ — —6 x 10*'" erg, 
which is comparable to the binding energy of the CD galaxy 
'I'CD ~ —(GMqj^/Rci))- It is however significantly smaller 
than the energy radiated from gas inside the cooling radius 
for the age of the cluster (e.g., for five Gyr this amounts 
to 8 X lO^'^erg). Nevertheless, this energy is easily avail- 
able from fast moving galaxies from beyond the very inner 
region — for example its value coincides with the kinetic en- 
ergy in a mass similar to that of the CD galaxy and moving 
at 2000 km s"'^, and that material can be supplied solely by 
the stellar mass of galaxies within 2rcooi — the region from 
which, according to our model, the bulk of energy input to 
the gas comes from. 

It is worth noting here that, under the circumstances 
just described, the mass of the gas within TcooI [which using 
equation (4) of Churazov et al. 2003, adjusted for h = 0.7 
with rcooi — 130 kpc amounts to 2.2 x IO^^Mq] is comparable 
to that of the galaxies in that region. Due to the tempera- 
ture drop in the gas, thermal motions can be significantly 
smaller than that of the galaxies. The energy of the gas can 
therefore be smaller than that of the galaxies by up to an 
order of magnitude. There is therefore sufficient energy in 
galaxies to heat the gas many times over, with the crucial 
parameter, which was the focus of this paper, being the rate 
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at which this is transferred. Furthermore, once the dark mat- 
ter in the center has been heated, it absorbs little additional 
energy from the galaxies, since dynamical friction does not 
act on fast particles, and the coupling with the dark matter 
at larger radii decreases rapidly with radius (see El-Zant et 
al. 2003 for further discussion). The energy lost to the dark 
matter is therefore less than that going into keeping the gas 
at constant temperature. The mass of the CD galaxy, as 
well as the spatial and velocity distribution of galaxies in 
the central region of the cluster will reflect the history of 
energy loss to both components. 

Several questions remain open. Prominent among these 
is the issue of thermal stability. From the condition that 
thermal stability requires a remarkably narrow range of 
heating rates, Bregman & David (1989) have argued against 
Miller's proposition that the Perseus cluster gas is heated 
by DF from galaxies. However the functional form of the 
heating rate used by these authors does not agree with 
later calculations that have been borne out by numerical 
studies — e.g., their postulated form has a negative energy 
transfer rate (~ ~l/V^) for highly subsonic velocities, re- 
sulting in gas cooling, which is incompatible with the posi- 
tive (even if small) heating found by Just et al. (1990) and 
Ostriker (1999) in that limit. It will be therefore necessary 
to repeat these calculations using adequate forms for the 
heat-loss function. These calculations would also address an- 
other crucial question, that concerning the precise manner 
in which the energy lost by the galaxies is distributed in 
the gas. We have adopted the simple approach where this 
energy is deposited isotropically and equally in logarithmic 
intervals. This should approximate how the energy is ini- 
tially deposited, at least for highly supersonic galaxies (and 
since, in the cooling region, the sound speed may be several 
times smaller than the velocity dispersion, this may not be 
a very bad assumption). Even then, however, the deposited 
energy may still be transported toward the center of the 
cluster by wave motions, as suggested, for example, by Bal- 
bus & Soker (1990), and redistributed. We note here that 
the claim of these authors that DF from galaxies is insuffi- 
cient to heat the cluster core is based upon outdated values 
for the electron density and an outdated (1/r) gas distribu- 
tion. With updated values for the electron density and gas 
distribution (e.g., the empirical formula of Churazov et al. 
2003), their conclusions are changed, as we have shown. We 
thus conclude that it is far from obvious that DF-heating of 
cluster gas is irrelevant. 



ACKNOWLEDGMENTS 

The authors would like to thank M. C. Begelman, E. 
Bertschinger, P. Goldreich, R. Narayan, E. Ostriker, N. 
Scoville and R. Sunyaev for helpful discussions and commu- 
nications. This work was supported at Caltech by NASA 
NAG5-9821, DoE DE-FG03-92-ER40701 and NSF grant 
AST 00-98301, and at the CfA by NAG5-10780 and NSF 
grant AST 0307433. 



REFERENCES 

Balbus S. A., Soker N., 1990, ApJ, 357, 353 
Bregman J. N., David L. P., 1988, ApJ, 326, 639 
Bregman J. E., David L. P., 1989, ApJ, 341, 49 
Briiggen M., Kaiser C. R., 2002, Nature, 418, 301 
Brunzendorf J., Meusinger H., 1999, A&AS, 139, 141 
Cho J., Lazarian A., Honein A., Knaepen B., Kassinos S., 

Moin S., 2003, ApJ, 589, L77 
Churazov E., Sunyaev R., Forman W., Bohringer H., 2002, 

MNRAS, 332, 729 
Churazov E., Forman W., Jones C, Bohringer H., 2003, 

ApJ, 590, 225 
Ciotti L., Ostriker J. P., 2001, ApJ, 551, 131 
Dokuchaev V. P., 1964, Sov. Astron., 8, 23 
El-Zant A. A., Hoffman Y., Primack J. P., Combes F., Shlos- 

man I., astro-ph/0309412 
Fabian A. C., 1994, ARAA, 32, 277 

Fabian A. C., Voigt L. M., Morris R. G., 2002, MNRAS, 
335, L71 

Gerhard O., Kronawitter A., 2001, AJ, 121, 1936 

Ghigna S., Moore B., Governato F., Lake G.,Quinn T., 

Stadel J., 2000, ApJ, 544, 616 
Just A., Deiss B. M., Kegel W. H., Bohringer H., Morgil G. 

E., 1990, ApJ, 400, 353 
Kaastra J. S. et al. , 2004, A&A, 413, 415 
Kaiser C. R., Binney, J., 2003, MNRAS, 338, 837 
Kim W.-T., Narayan R., 2003, ApJ, 596, 889 
Kim W.-T., Narayan R., 2003, ApJ, 596, L139 
Miller L., 1986, MNRAS, 220, 713 
Narayan R., Medvedev M. V., 2001, ApJ, 562, L129 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 

493 

Ostriker E., 1999, ApJ, 513, 252 

Peres C. B., Fabian A. C, Edge A. C., Allen S. W., John- 
stone R. M., White D. A., 1998, MNRAS, 298, 416 
Peterson J. R. et al. , 2001, A&A, 365, L104 
Peterson J. R. et al. , 2003, A&A, 590, 207 
Raphaeh Y., Salpeter E. E., 1980, ApJ, 240, 20 
Ruderman M. A., Spiegel E. A., 1971, ApJ, 151, 679 
Salucci P, Ashman K. M., Persic M., 1991, ApJ, 379, 89 
Sanchez-Salcedo F. J., Brandenburg A., 1999, ApJ, 522, L35 
Sanchez-Salcedo F. J., Brandenburg A., 2001, MNRAS, 322, 
67 

Takamiya T., Sofue Y., 2000, ApJ, 534, 670 

Trujillo I., Burkert A., Bell E. F., 2004, ApJ, 600, L39 

Tucker W. H., Rosner, R., 1983, ApJ, 267, 547 

Voigt L. M., Schmidt R. W., Fabian A. C., Allen S. W., 

Johnstone R. M.n 2002, MNRAS, 335, L7 
Voigt L., Fabian A., 2004, MNRAS, 347, 1130 
Zakamska N. L., Narayan R., 2003, ApJ, 582, 162 



© 0000 RAS, MNRAS 000, 000-000 



